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Abstract 

We show that long-term memory effects, present in the chaotic disper- 
sion process generated by a meandering jet model, can be nonetheless 
taken into account by a first order Markov process, provided that the 
states of the phase space partition, chosen to describe the system, be 
appropriately defined. 

1 Introduction 

Geophysical processes often exhibit an irregular behavior that is due to de- 
terministic chaos and/or to the presence of many relevant degrees of freedom. 
Striking examples are climate dynamics (Praedrich 1988, Nicolis et al. 1997) 
and transport processes in the atmosphere and in the oceans where many 
different characteristic temporal scales are involved (Samelson 1992, Yang 
1996). After Lorenz's seminal work (Lorenz 1963), it is now well recognized 
that even deterministic systems with as few as three degrees of freedom 
can have complex temporal evolution similar to genuine stochastic processes 
(Ott 1993, Beck and Schlogl 1993). For example, the large-scale and long- 
time dispersion characteristics of drifters deterministically advected by a non- 
stochastic velocity field, giving rise to Lagrangian chaos, can, sometimes, be 
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well described by random walk models (Ottino 1989, Crisanti et al. 1991). 
On the other hand, a simple probabilistic description, e.g., one in terms of 
autoregressive models and Markov chains (Vautard et al. 1990, Nicolis et al. 
1997), is often not able to catch the richness of most geophysical phenomena. 

In a recent paper (Cencini et al. 1999), mixing and transport in a mean- 
dering jet were investigated and compared with a Markov chain. It was found 
that, due to the presence of very different characteristic time scales, this sim- 
ple approach was not fully satisfactory. In fact, the ubiquitous presence of 
a hierarchy of time scales forces one to develop nonstandard techniques for 
data analysis and for their modeling. 

The main aim of this Note is to show how a suitable Markov chain with 
a small number of states can reproduce in a proper way the main statistical 
features of systems similar to the one that is considered below. The basic 
point is an appropriate definition of the states of the system, such that long- 
time memory effects are straightforwardly taken into account in a Markov 
chain approximation of the dynamics. We expect this kind of analysis to be 
of interest in the more general context of systems with multiple time scales. 
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2 Transport in a meandering jet 

The test case under consideration consists of a 2D meandering jet, formerly 
introduced as a model for studying Lagrangian transport across the Gulf 
Stream (Bower 1991, Samelson 1992). In this model, the steady flow is 
described by the following stream function: 

, I , V — B cos (kx) I , , 

y) = - tanh \ fc-J=_4=^ \ + cy (1) 

where x and y are the fluid particle zonal and meridional coordinates, re- 
spectively, k is the wavenumber of the flow along the x-direction, B controls 
the amplitude of the meanders and c is the retrograde, counter stream, ve- 
locity component far from the jet's core which is centered on y = 0. The 
advective time, Tadv, is defined as the time needed for a particle, traveling 
close to the jet core, to move along the x direction a distance equal to the 
wavelength of the flow A = 2Tr/k. For the typical values of the parameters 
used in our calculations, the value of this characteristic time is Tadv — 47r. 
The recirculation time. Tree, is defined as the period of a trajectory moving 
on a closed orbit next to the boundary of a gyre. This time is Tree 10' T^dv • 
The 2D incompressible velocity field {u,v) and the stream function ^(x,y) 
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are related by, 



dy' dx' 



The flow is invariant under north-south inversion {y direction) and is periodic 
along the x-axis with period A. In the absence of perturbations, passive 
particles follow the streamlines at constant \E'. The steady flow corresponding 
to this stream function consists of open streamlines, to be identified as "the 
jet", and closed streamlines, to be called "the gyres", see Figure |I[ Non- 
trivial behavior appears as soon as a perturbation is added to the above- 
described flow. In particular, cross-stream transport, i.e., from gyre to gyre, 
occurs either by the action of a random forcing or due to non-stochastic, time 
dependent terms inserted in the stream function which lead to the onset of 
Lagrangian chaos. In the stochastic model the position of a passive particle 
evolves according to 

^ = u{x,y) +wi{t), ^ =v{x,y) +W2{t) (3) 

where the stochastic flelds Wi(t) have red-noise correlations, i.e., 

< Wi{t)wj{t') >= 2a25i,e-'*-*'l/^ t,j = 1, 2, (4) 

with 0"^ the variance of the stochastic flelds and r their correlation time. 
The deterministic, chaotic model is obtained by making the amplitude of the 



meanders, B in equation (p, time dependent, 

B =^ Bit) = Bo + e cos(cjt + (f)) (5) 

where Bq is the mean meander amphtude and e and uj are appropriately 
chosen so as to generate large-scale mixing, as in Cencini et al. (1999), i.e., 
leading to an overlap of resonances (Chirikov 1979); (p is an arbitrary phase. 

The symbolic dynamics (Beck and Schlogl 1993) of the noisy model can be 
satisfactorily approximated by a simple minded, first order Markov process. 
In the sequel we show that, in order to describe the chaotic djTiamics via a 
first order Markov chain, it is necessary to introduce a more suitable definition 
of "state" which take into account the non trivial dynamical behavior of the 
system. This is made clear by comparing the symbolic dynamics generated 
by two different partitions of the phase space, both partitions having four 
cells. 

The simple minded partition, call it 11, is defined by the southern gyre 
(symbol 1), the southern half jet (symbol 2), the northern half jet (symbol 
3) and the northern gyre (symbol 4), see Figure |^. The partition takes into 
account the north-south symmetry as well as the periodicity of the system. 
The position of the advected particles is sampled millions of times at a time 
interval T = 27t/ uj. At each sampling time, one associates one of the symbols 
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according to the partition cell in which the particle is located at that time. 
In this way, a symbolic sequence is generated from each trajectory. Once this 
has been done, one computes the transition matrix W, giving the state-to- 
state transition probabilities in the time interval T. 

While the precise nature of the perturbation, i.e., stochastic or purely 
deterministic, does not affect the gross qualitative features of the dispersion 
process, it does lead to very different correlation functions. In the determinis- 
tic, chaotic case, the transport process shows persistent, nontrivial long term 
correlations which are suppressed when stochastic noise is present. These 
differences are clearly seen in Figure |^. The Wij matrix element is the prob- 
ability to observe state j at time {n + 1)T knowing that the state i occurs 
at time nT (where n is integer). In the present case, using partition 11, this 
analysis leads to. 



W 



( \ 
0.851 0.147 0.002 0.000 



0.247 0.535 0.215 0.003 



0.003 0.215 0.535 0.247 



y 0.000 0.002 0.147 0.851 j 



(6) 



The symmetries in the transition matrix W reflect the dynamical and sta- 



tistical equivalence of the northern and southern halves of the system. For 

future reference, we also quote the corresponding steady probability vector 

P = (Pi, -P2, -P3, Pa) which satisfies the equation 
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1=1 

One finds: in the gyres, P\ = P^ = 0.315, and in the jet, P2 = -P3 = 0.185. 

In contraposition to the stochastic model defined by Eqs. (^) and (^, 
the deterministic evolution with (^) is characterized by long-term correla- 
tions which are poorly reproduced by the Markov chain generated from W. 
Essentially, the matrix elements and overestimate the coupling be- 
tween the two halves of the systems. This failure of the simple probabilistic 
model based on W is illustrated in Figure ^ where the actual autocorrelation 
functions of the states 1 and 2 are compared with the corresponding corre- 
lation functions generated by the first order Markovian process described by 
W. 

As one can see, in the actual dynamics there are basically two character- 
istic times, Tf and r^, that correspond to the inverse of the two exponential 
decay rates at short and at long times, respectively. As order of magnitude, 
the shorter decay time is tj ~ T, and the longer one is ~ 10^ ■ T. 

A way of estimating a lower bound for the order of the Markov chain which 
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would be necessary in order to adequately reproduce the system's statistical 
features, see, e.g., (Khinchin 1957), is based on the block entropies Hn which 
are defined by 

Cn 

where P{Cn) is the probability of observing a sequence C„ = {it+i,it+2, , it+n) 

of n successive symbols generated by the dynamical system, with e 
{1,...,4} the partition cell visited at the time t = mT. From the block 
entropies if„ one computes next the quantities hn — Hn — Hn-i which rep- 
resent the average additional information needed to specify the n-th symbol 
in given the sequence C„_i. The limit of /i„ for n — > oo gives the Shannon 

entropy hs of the infinite sequence ii,i2, For a Markov process of order 

i/, one has hn = hs for all n > u + 1, (Khinchin 1957). In agreement with 
what was found in the paper by Cencini et al. (1999), at least in order to 
reproduce the entropic properties, one would need a Markov process of very 
high order. Recall that this procedure gives only a lower bound to the order 
of the Markov chain. It is clear that the strategy consisting of simply in- 
creasing the order of the Markov process is not practical at all. Surprisingly, 
we were able to give a non-obvious choice of the partition that allows us to 
work with a first order Markov chain where now the time memory effects are 
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satisfactorily described. This is explained in the following. The new parti- 
tion, denoted by 11*, is a 4-state partition similar to the previous 11 but the 
jet is partitioned in a different way: the jet corresponds to symbol 2 if the 
last previously visited cell is the southern gyre 1 or the cell 2 itself; if the last 
previously visited cell is the northern gyre 4 or the cell 3 itself, then the jet 
corresponds to symbol 3, see Figure |^. Accordingly, transitions between the 
two jet states 2 and 3 are not possible. This partition of the jet preserves 
part of the long time memory of the gyre-to-current transitions, a feature of 
the deterministic chaotic model which is impossible to reproduce with the 
previous partition 11. In other words, this partition introduces some memory 
even if formally the model remains of first order. The time signal generated 
by the 11* Markov process is compared to the actual time signal in Figure ^. 
The memory effects are visible through the existence of "blocks" describing 
the fast back and forth gyre-to-current transitions in one half of the system, 
before moving to the other half. The new transition matrix W* is found to 
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be, 



W* 



( \ 
0.850 0.150 0.000 0.000 



0.243 0.749 0.000 0.008 



0.008 0.000 0.749 0.243 



y 0.000 0.000 0.150 0.850 y 



(8) 



The corresponding time-independent probability vector P* is equal to the 
one found for the 11 partition, P* = P, as it should be. Notice that, because 
the transitions between states 1 and 2 (3 and 4) are much more frequent than 
those from state 3 to 1 (from 2 to 4), the matrix W* is not far from being 
block-diagonal. 

The improvements introduced by the 11* partition show up in the cor- 
relation functions it generates, as can be seen in Figure |^, and in the fast 
convergence of the hn entropies shown in Figure H. Moreover, with the 11* 



partition the two characteristic times are well estimated. 



3 Conclusions 

We have shown that the long-term memory effects which are present in the 
chaotic dispersion processes generated by the flow defined by Eqs. (|l]) and 
d^) can be reproduced by a Markov chain, provided that the partition of the 
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phase space is performed in a special way. It is important to remark that 
first order in a Markov process generated with the 11* partition corresponds 
to a very large order, say ~ 0{ts/T), in a Markov process generated from the 
n partition. We hope that the ideas presented in this Note will be of use in 
a more general context and we plan to apply the above discussed technique 
to the analysis of other geophysical phenomena. Needless to say, when more 
than two characteristic times are relevant, partitions with a larger number 
of cells will be required. 
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FIGURE CAPTIONS 
FIGURE |1]. Schematic diagram of the circulation in the meandering jet 
model. Arrows indicate the direction of particle motion: "westerly" ad- 
vection in the jet current, clockwise and anti-clockwise recirculation in the 
southern and northern gyres, respectively. Spatial coordinates X and Y are 
given in unit of the wavelength A = 27r/A; and Bq is set equal to 0.16A, see 
eq. (0). 

FIGURE ||. n partition of the steady meandering jet: 1) south gyre; 2) south 
half-jet; 3) north half-jet; 4) north gyre. 

FIGURE^. Comparative plot of the two typical symbolic sequences obtained 
from: a) the stochastic model and b) the chaotic model, with the parameters 
in eq. being e = Bq/4 and u ~ 2TT/Tadv The parameters of the red noise 
in eq. (D are = 6 ■ lO~^{X/Tadv f and r = T/4. 

FIGURE |[ Autocorrelation functions for the states 1 and 2 of the 11 par- 
tition, see Figure ^: actual autocorrelations from the chaotic model (full 
lines) and those generated by the transition-rate matrix W (dotted hues). 
Parameters of the deterministic perturbation: as in Figure |^. 
FIGURE |[ n* partition of the steady meandering jet: 1) south gyre; 2) and 
3) jet; 4) north gyre. The central jet corresponds to state 2 if the last visited 
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gyre is state 1, and to state 3 if the last visited gyre is state 4. 
FIGURE ^ Comparative plot of: a) the chaotic model symbolic dynamics 
using the partition 11* and b) the symbol sequence simulated on the basis of 
the transition-rate matrix W*. 

FIGURE ^ Autocorrelation functions for the states 1 and 2 of the 11* parti- 
tion, see Figure ^: actual autocorrelations from the chaotic model (full lines) 
and from the Markovian approximation generated with W* (dotted lines). 
FIGURE |[ Differential Block Entropy hn = Hn+i — Hn, as a function of the 
sequence-length n, for the 11 (*) and the 11* (□) partitions. For each length 
n approximately ~ 10'' sequences were used. Notice the very different rate of 
convergence of hn to its limit value (Shannon entropy) /is in the two cases. 
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Figure 1: 
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Figure 2: 
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